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We apply a variational method devised for the nuclear many-body problem to the 1-dimensional 
Hubbard-model with nearest neighbor hopping and periodic boundary conditions. The test wave 
function consist for each state out of a single Hartree-Fock determinant mixing all the sites (or 
momenta) as well as the spin-projections of the electrons. Total spin and linear momentum are 
restored by projection methods before the variation. It is demonstrated that this approach repro- 
duces the results of exact diagonalisations for half-filled N = 12 and N = 14 lattices not only for 
the energies and occupation numbers of the ground but also of the lowest excited states rather well. 
Furthermore, a system of 10 electrons in a N = 12 lattice is investigated and, finally, a N = 30 
lattice is studied. In addition to energies and occupation numbers we present the spectral functions 
computed with the help of the symmetry-projected wave functions, too. Also here nice agreement 
with the exact results (where available) is reached. 
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I. INTRODUCTION 

Consider a finite number of TV identical Fermions in a model space defined by some finite number M of suitably 
chosen orthonormal single Fermion basis states. Assume furthermore, that the effective Hamiltonian of these Fermions 
appropriate for this model space is known. Then, at least in principle, this many-Fermion problem is exactly solvable. 
One only has to distribute the N Fermions over the M orbitals according to the Pauli principle, i.e., construct all 
possible Slater-determinants, and then to diagonalize the known Hamiltonian in the resulting configuration space of 
dimension (^). 

The nuclear shell-model (see Q for a recent review and references therein) is one realization of this scheme. In order 
to describe ground-state and low energy excited states of a nucleus Z v valence protons and N v valence neutrons are 
distributed over M p and M n single-particle basis states and a suitably chosen effective Hamiltonian is then diagonalized 
in the resulting configuration space. The dimension of this configuration space which is (^ p ) • ( N n ) , can, however, 
be reduced considerably, if use is made of the symmetries of the effective Hamiltonian. So, e.g., the effective nuclear 
many-body Hamiltonian is a scalar in normal space and thus does neither mix states with different total angular 
momentum nor with different z-projections of the latter. Furthermore, neglecting weak interactions, parity becomes 
a good quantum number, too. Consequently, the above Slater-determinants can always be coupled to configurations 
with definite parity and angular momentum quantum numbers and the Hamiltonian can then be diagonalized in the 
much smaller configuration spaces corresponding to such a symmetry representation. However, even if use is made of 
all the symmetries, in general the dimensions are far too large to be numerically tractable. Thus in most cases one 
has to rely on approximate methods, which truncate the complete shell-model expansion of the wave functions to a 
numerically feasible number of configurations. How this can be done without loosing the essential degrees of freedom 
being relevant for any particular state under consideration is the central question of nuclear structure physics. 
The Hubbard-model, a schematic model developed to describe some basic features of solid state physics, is another 
example for the scheme outlined above. In this case N e electrons are distributed over the N sites of a lattice and 
the corresponding model-Hamiltonian is then diagonalized within the resulting configuration space (see, e.g. Q and 
references therein). Again, the dimension of this configuration space increases drastically with the number of lattice 
points N and again, even if all the available symmetries are used, the dimensions are in general far too large to allow 
for an exact solution. Thus again approximate methods are asked for and rather successful ones have been developed 
within the past Q. The aim of the present study is to adopt an approximation scheme, which has successfully been 
applied in nuclear physics, and apply it to the Hubbard model to test its efficiency for this system. 
One of the most successful truncation approaches to the nuclear many-body problem has been addressed in a recent 
review article jjj. It works with general Hartree-Fock-Bogoliubov configurations, breaking all the symmetries required 
by the effective nuclear many-body Hamiltonian, as basic building blocks. These symmetries are restored with the 
help of projection techniques and the underlying HFB-transformations as well as the configuration-mixing are then 
determined by variational calculations. In quantum chemistry such kind of approaches are called unrestricted self- 
consistent field approximations with variation after projection and it is known that electronic correlations can be 
accounted for this way by breaking symmetries 
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Here we want to address the question how well the symmetry projected variational approach will perform when 
applied to the Hubbard model. For simplicity we shall restrict ourselves to the 1-dimensional Hubbard model and use 
only general Hartree-Fock instead of Hartree-Fock-Bogoliubov configurations in the symmetry projected variational 
approach. As a test we will determine energies and single-particle spectral functions and compare the results to exact 
results for cases, for which the exact diagonalization can be done. We will also demonstrate the feasibility of the 
symmetry projected variational approach for cases, which cannot be solved by brute force diagonalization. 
We present the relevant formalism in sec. [H] We start with a short sketch of the 1-dimensional Hubbard model 
with nearest neighbor hopping and periodic boundary conditions in sec. Ill Al In sec. Ill Bl we shall then derive the 
variational equations for the Hartree-Fock approach with spin- and linear momentum-projection before the variation 
for the ground as well as for the excited states of the considered system. How to calculate the corresponding (hole- a 
well as particle-) spectral functions is discussed in sec. Ill CI Sec. lIIII contains some general properties of the Hubbard- 
model and presents then the results obtained within the variational approach for half-filled TV = 12-, TV = 14- and 
TV = 30-lattices. Where possible these results are compared with those of complete diagonalisations. Furthermore, 
some results for a 10 electron system in a TV = 12 lattice are presented. Finally, some conclusions and an outlook are 
given in sec. HVI 



II. THEORY 



A. The 1-dimensional Hubbard model 

In its simplest version Q the effective Hamiltonian of the 1-dimensional Hubbard-model has the form 

JV +1/2 ( \ N 




cjaCj+la > + UJ2 cjTSV^T ■ (1) 
J j=l 

Here the operator c| CT creates from the particle vacuum |0) an electron with spin-projection a = ±1/2 along an 
arbitrary chosen quantization axis on the site j (j = 1, . . . , TV), the corresponding annihilator Cj a destroys such an 
electron. Obviously, these operators fulfill the standard anti-commutation relations for Fermion-operators 

{4<t' £ J' CT '} = d ]a' c 3'<y' + Cfo-'cJ* = 5 j,j' S <r,<r> ■ ( 2 ) 

The Hamiltonian Q simulates a system of electrons in a periodic potential, in which each of the TV potential wells is 
supposed to have only a single bound electron state, which can be occupied by at most two electrons with opposite 
spin-directions. In case that two electrons are occupying the same state, they feel the repulsive Coulomb-interaction 
(U > 0). Furthermore, each electron may tunnel to the neighboring well (if the corresponding site is not already 
filled by an electron with the same spin-projection). This so-called "nearest-neighbor-hopping" is described by the 
"hopping parameter" t > 0. As usual, all energies are measured in units of this parameter and TV is supposed to be 
an even integer. One furthermore assumes periodic boundary conditions, i.e. the sites TV + 1 and 1 are identical. 
Thus the system lives on a circle of length L = TV A. The spacing A is set to unity in the following. Let us for the 
moment consider a "half-filled" grid, i.e. N e = TV electrons on the TV sites. It is obvious, that for interaction strength 
U = these electrons form a noninteracting "Fermi-gas" in a finite, 1-dimensional box, while for very large interaction 
strength (U — * oo) an anti-ferromagnetic ground state with total spin S = is to be expected Q. 
We apply now the Fourier transformation 



fit 

Qf(T 



7= J2 CX P{-*^ j} c\ a (3) 

on the basis orbits (A = 1). This yields a set of TV single electron states in momentum space with momenta 

k a = ^a a = -TV/2 + 1,. ..,N/2, (4) 
in terms of which the Hamiltonian gets the form 

1/2 

a cr=-l/2 ^ ' a, 0,7, 5 
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where the generalized Kronecker-symbol is one, if a + (3 — 7 — S is either or ±N (because of the periodic boundary 
conditions), and zero else. The linear momentum quantum numbers a, /?, 7, S run all over —N/2 + 1, . . . N/2. Because 
of the "cosine-dispersion" in the one-body term, the single particle spectrum contains one state with energy — 2t (for 
a = 0), one state with energy +2t (for a = N/2) and N — 2 two-fold degenerate states with energies cos([2ir / N]a) 
(fora = ±l,... ) ±(jV/2-l)). 

Obviously, the Hamiltonian (QJ or (0 conserves the total number of electrons. Furthermore, it is easy to show, that 

» 2 

it commutes with the square of the total spin operator 5 as well as with its 3-component S z and is thus a scalar 
in spin-space. Hence its eigenstates can be classified according to the corresponding spin-quantum number 5 and 
are degenerate for all the 25+1 values of its z-projection £ = —S,...,3. For an even number of electrons only 
integer values for the total spin can occur. Since each of these spin-values has aE = component, it is sufficient to 
diagonalize or © in the space of the S = configurations in order to obtain all the eigenstates. In contrast, for 
an odd number of electrons 5 is a half-integer, so that in this case all the configurations with £ = 1/2 have to be 
used as basis states. The lowest eigenvalues and eigenvectors for such a matrix with large dimension can be obtained 
by applying efficient algorithms like e.g. the Lanczos-method 0. For half-filling (N e = N even) the total number of 
S = configurations to be treated is 

n(N/2, S = 0) = ( A N /n ) , (6) 



\N/2 

and thus increases drastically with the number N of available sites (for TV = 16, e.g., eq. © yields already 165 636 
900 configurations). Hence the applicability of this approach is rather limited. The situation becomes slightly better, 

if instead of all £ = Slater-determinants only those configurations with a definite total spin-value 5 are included. 
Again for half-filling, one obtains here dimensions of 

"(»m s U- s ) ! -U-Vi)' < 7 » 

which for N = 16 still amount to 34 763 300 5 = and even 66 745 536 5=1 configurations, respectively. In 
addition, such a procedure requires the coupling of the simple Slater-determinant to configurations with good total 
spin and thus complicates the calculation of the various matrix elements considerably. 

There is, however, a further symmetry of the Hamiltonian, which can be used to reduce the dimension of the configu- 
ration spaces. As can be shown easily, the operator commutes (modulo +2ir) with the Hamiltonian of total linear 
momentum 



P = Y,\^ a \ £ ^aa , (8) 



a, a 



and thus the eigenstates can be classified in addition by a momentum quantum number £, too. The corresponding 
total momentum is = (27r/TV)£. The total linear momentum quantum number can assume all values given in eq. <@} 
or, equivalently, all integer values (£) in between and (N — 1), where use has been made of the periodic boundary 
conditions. In addition, it is immediately seen, that the Hamiltonian (01 remains unchanged, if the signs of all the linear 
momentum quantum numbers are flipped (because of the periodic boundary conditions —N/2 — > —N/2 + N = N/2). 
Thus the states with £ = 1 are degenerate with those with £ = JV — 1, those with £ = 2 with those with £ = N — 2, 
etc., while the states with £ = TV/2 (again, since —N/2 — > —N/2 + N = N/2) occur only once. 

For half-filling, the total number of configurations with a given total spin 5 and given linear momentum quantum 
number £ cannot be given in a closed form but it can be easily calculated numerically. It turns out that 

n(N/2, 5£) ~ n(N/2, S)/N (9) 

is a rather good first guess for all possible £-values, if N is sufficiently large. For N = 16 and 5 = 0, e.g., the dimensions 
for the various configuration spaces are 2 172 400 (for £= 1,3, 5, 7, 9, 11, 13, 15), 2 173 008 (for £ = 2, 6, 10, 14), 2 
173 016 (for £ = 4, 12) and 2 173 018 (for £ = 0, 8), respectively. 

However, even if spin and linear momentum are used as "good" quantum numbers, the dimension of the resulting 
configuration spaces still increases drastically with the number of sites. For TV = 18 and 5 = eq. © yields already 
2.5 * 10 7 , and for TV = 30 and 5 = even 9.7 * 10 13 configurations for each £ have to be treated. Thus, at least for 
half-filling and TV > 18, an exact solution is almost impossible even on modern computers. This statement remains 
valid even if we furthermore consider the symmetry of the Hubbard Hamiltonian under charge rotation and the 
corresponding conserved quantum numbers. 
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Therefore one has to rely on approximate methods, which truncate the complete expansion of the wave functions to a 
numerically feasible number of configurations. How this can be done without loosing the essential degrees of freedom 
relevant for the particular states under consideration is the central question not only for the Hubbard-model but for 
any other finite many-body problem as well. 

One of the most favored methods for such a truncation is provided by Greens-function Monte-Carlo calculations, in 
which the relevant configurations are selected according to their statistical weight. This method has been applied 
rather successfully not only to the Hubbard-model @ but also to the nuclear many-body problem 0. It has, however, 
the limitation that only expectation values within the ground state or in some statistical assembly can be obtained. For 
Fermion-systems the method suffers in addition from the well-known "sign-problem". In nuclear physics both these 
shortcomings could be overcome in recent years by the so-called "Quantum-Monte-Carlo-Diagonalization" (QMCD) 
method Q, in which again the relevant configurations are selected by stochastic methods, however, the sign-problem 
can be bypassed and ground as well as excited states can be treated on equal footing. 

However, stochastic selection is not the only promising method for truncation in finite many-body problems. Another 
possibility is provided by variational methods. Here one uses the most general Slater-determinants (or even gener- 
alized Slater-determinants), which can be constructed within the considered single particle space, as basic building 
blocks. Unfortunately, in general these configurations do break all the symmetries required by the chosen many-body 
Hamiltonian and hence cannot be considered as physical states but only as an approximation introduced to account 
for as much as possible of the correlations in as few as possible configurations. The required symmetries, however, can 
be restored with the help of projection techniques, and the resulting symmetry-projected configurations can then be 
used as test wave functions in chains of variational calculations in order to determine the underlying single particle 
transformations as well as the configurations mixing. Such symmetry-projected variational calculations on the basis 
of general Hartree-Fock-Bogoliubov configurations have been applied very successfully to the nuclear many-body 
problem within the last two decades (see Q for a recent review and references therein), and it could be proven that 
they work equally well as alternative approaches like, e.g., the QMCD method. 

In the present work we want to demonstrate that these symmetry-projected variational methods are not only useful for 
the nuclear many-body problem but can also be applied to other finite many-Fermion problems like the 1-dimensional 
Hubbard-model described above. 



B. Symmetry projected Hartree— Fock 

For this purpose we start by introducing "quasi-particles" of the type 

n/2 +1/2 

" b * = E E Ka, a cL> (10) 

a=-iV/2+l ct=-1/2 

where D is a linear (2N x 27V) transformation, which has to be unitary 

D*D = DD^ = 1 2N , (11) 

in order to conserve the Fermion anticommutation relations for the quasi-particle creators ijlOl) and the corresponding 
annihilators. Eqs. iflOll and l(TT|l describe a general Hartree-Fock (HF) type transformation. It should be stressed 
here, that eq. lllOt could still be generalized by including linear combinations of the annihilators on the right hand 
side. This would result in a so-called Hartree-Fock-Bogoliubov (HFB) transformation as used, e.g., in the approaches 
reviewed in ref. 0. In the present paper, however, we shall restrict ourselves to the simpler HF-transformations only. 
Note, that nevertheless the transformation lllOt mixes all the linear momentum quantum numbers as well as the 
spin-projections of the basis states <EJ). 

In the usual Hartree-Fock approach one searches then for the optimal one-determinant representation of the N e - 
electron ground state 

\°) = In s ui°>' ( i2 ) 

in which the energetically lowest N e states b^ (h = 1, . . . , N e ) of the form IjlOJI are occupied and the remaining 2N—N e 
states bj, (p = N e + 1, . . . , 2N) are empty. In the following the notation h, h', . . . is always used for the occupied 
states, while the unoccupied ones will be denoted by p, p',. . ., respectively. 

Obviously, the determinant l(T2|l conserves the total number of electrons (N e ) but is neither an eigenstate of the 
square of the total spin operator S nor of its 3-component S z nor of the total linear momentum operator {BJ . These 
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symmetries have therefore to be restored with the help of projection techniques. For the spin-quantum numbers this 
can be achieved via Villars' famous projection operator 



P lv> = J dQD% x ,(Q)Rs{n), (13) 

where Rs(ty is the rotation operator in spin-space, the Wigner-function Z?| E , (fi) its representation in spin- 
eigenstates and the integral has to be taken over the full range of the three Euler-angles. Because of the non-abelian 
nature of the rotation group, il 1 31) is not a true projector in the strict mathematical sense. In order to achieve 
independence of the choice of direction for the "intrinsic" quantization axis one is forced to use the linear combination 

+s 

\D;N e SX) = P lv\D)fc> (14) 

E'=-S 

as test wave function in the variation with the mixing coefficients / to be treated as additional variational variables. 
In a similar way the total linear momentum can be restored. This is done via the operator 

tf(0^E^{i(*-f*)i}. (15) 

which projects the determinant H12I1 on the component with linear momentum k = (2-k/N)^. The operator \1q \ is 
the finite, 1-dimensional limit of the general operator restoring Galilean invariance discussed, e.g., in ref. Q. In 
contrast to nuclear systems where the Galilean principle of relativity imposes k = 0, lattice systems allow solutions 
with k > 0: the Hamiltonian QJ or © has to be considered on the (infinitely heavy) background of the ions providing 
the periodic potential. This background can absorb any change of linear momentum of the electrons easily so that 
Galilean invariance for the total system is always ensured. It should be stressed furthermore, that the spin-projection 
and the linear momentum projection have to be performed before the variation. Then the correct moment of inertia 
and mass is restored [lfij . 

Using 115t in addition to llllll we obtain the projected determinant 

+s 

\D;N e ^SE}= E Pi,x>C(0\D) fc> (16) 

as ansatz for our test wave function. The corresponding energy functional 

r _ (D; N^S(S)\H\D; JV e £S(£)) 



(D; N e £SC£)\D; JV f5(E)) 



(17) 



where the spin-projection £ has been put in parentheses since the total energy does not depend on this quantum num- 
ber, has now to be minimized with respect to the mixing coefficients / as well as to the underlying HF-transformation 
D. 

Variation with respect to the /'s yields the generalized eigenvalue problem 

(H-EN)f = 0, (18) 

with the constraint 

pNf = 1 2S+1 (19) 
ensuring the orthonormality of the solutions, and the (25"+ 1) x (2S+ l)-dimensional square matrices N and H given 

by 

7V E . S , (D\P£^,C(0\D) (20) 

and 

H^, = (D\HPi )S ,C(0\D) , (21) 

respectively. The matrices N and H are hermitian and the overlap matrix ll2f))l is furthermore positive definite. In 
eqs. l(2fljl and l(2"T)) the obvious quantum numbers S and £ have been suppressed. In the following only the energetically 
lowest solution of eq. I|18H is kept. 



The minimization of the energy functional (I17|l with respect to arbitrary variations of the underlying HF- trans- 
formation D is more involved. This transformation has to be unitary and thus not all of the 2N x 2N matrix 
elements of D are linearly independent. Nevertheless, an unconstrained minimization of the functional itTTjl can still 
be performed, if one parameterizes the underlying HF-transformation via Thouless' theorem ^l|]> which states that 
any HF-determinant \D d ) can be represented in terms of the creation and annihilation operators of some reference 
determinant \Dq) via 

\D d ) = c(d)cxp d ph bl(D )b h (D ) j | Do) , (22) 

provided that the two determinants are non-orthogonal, since 

c(d) = (D \D d ). (23) 

The creation operators belonging to the HF-determinant \D d ) are then related to those of the reference determinant 
I -Do) via 

b{(Dd) = Y\L-\ h ' U{,(D ) + £ dp'h>bl,(D )\ (24) 



for the occupied and 



%(p d ) = J2i M '% P ' fe(A>) - £ d;, h ,b{,(D )\ 

V' V h> ) 



(25) 



for the unoccupied states, respectively. They are now given in terms of the (2N — N e ) ■ N e linear independent variables 
d p h- The lower triangular (N e x N e ) matrix L in J 2411 is defined by the expression 



LJV e 



d T d* = LL\ (26) 



while the lower triangular ((2N — N e ) x (2N — N e )) matrix M out of II25II can be obtained from the solution of the 
equation 

l 2 w-w e + d*d T = MM 1 . (27) 
Both, eqs. Il2(it and l|27ll . are usual Cholesky decompositions. 

The variational equations resulting from the minimization of the energy functional ifTTIl with respect to the HF- 
transformation now assume the form 



dE 



M~ l] G l- 1 



ee 0, (28) 

ph 



dd ph 

where the ((27V — N e ) x N e ) matrix G is defined by 

G ph = £ /! (D\(H - E1)P£^C(0 %(D) b h (D)\D) . (29) 



E. S' 



Once one has reached the solution not only the "global" gradient vector i|28|l but also the "local" one lt29jl does vanish 
identically. This vanishing of the gradient vector il29)l is a sort of "generalized Brillouin theorem" [Tl| : it describes 
the stability of the symmetry-projected HF-solution with respect to arbitrary symmetry-projected one particle-one 
hole excitations. 

For any given 5* and £, the simultaneous solution of the set of eqs. (tl8)) . (tl9)l and <l28)) yields the optimal representation 
of the energetically lowest state by one single symmetry-projected HF-type configuration. This corresponds to the 
so-called VAMPIR (Variation After Mean-field Projection In Realistic model spaces) approach out of ref. Q, though 
restricted here to a HF-type instead of HFB-type transformations. Since D (and hence d) as well as / are essentially 
complex matrices this solution results from the minimization of a function of m = 2 ■ (2N — N e ) ■ N e - real variables, 
if even N e and S = is considered, while in general for S > and arbitrar y N e , (AS + 1) + m real variables have to 
be treated. For this minimization a quasi-Newton procedure (see, e.g., [12lll3|l is used. 
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Excited states with the same quantum numbers S and £ can be treated in a similar fashion, if one ensures the 
orthogonality with respect to the solutions already obtained. For the first excited state (n = 2) this can be achieved 
with the help of the projection operator 



ra-1 



f„_! = J2 NetSEXD*; N e £SE\ , (30) 



which projects on the lowest solution (D 1 , f 1 ) of the form (tTfi)) . For the variational calculation then instead of l(T6jl 
we use its complement 

+s 

\D n ;N e £SX} ee (1-f) Pi,vC(t)\D n ) & (31) 

E'=-S 

as a test wave function for the variational calculation. The procedure ffiljl . lj3lt can be repeated for the second excited 
state (n = 3), etc., up to the n lowest (orthonormal) states. Finally, the residual interaction between these n states 
is diagonalized. This corresponds closely to the EXCITED VAMPIR approach from ref. Q. 

Furthermore, we would like to stress, that if a description by one single symmetry-projected determinant for each state 
is not sufficient, further correlations can easily be accounted for by successive variational calculations for additional 
determinants as done in the FED (from FEw Determinants) EXCITED VAMPIR approach out of ref. jjj. This, 
however, has not been done in the present work since one determinant gives already a good accuracy, as we will see 
below. 

Left to be computed are now the symmetry projected matrices lf20j) and (l2l)) generalized to two different determinants 
\D l ) and \D k ) on both sides, because of the eventual use of ll3"Tl) instead of ltl6|l as test wave function and of more 
than one determinant for the description of the N e ± 1-electron systems as discussed in sec. Ill CI Furthermore, the 
corresponding gradient vectors occurring in eq. I|29ll have to be calculated. One obtains successively 

(D i \PixAQ\D h ) = ^X>p(~^j^i^ / df2£>C s ,(QK fc (fi, j), 



■ cx p | 


f 2tt 


j=i 






Yl exp 


f 2tt 




ik (n, j)n 


"(Si, j 




(33) 



and, 



■n; k h (n,3)n* k (n :J ) : (34) 
(n l \HPl s ,C(Obl(D k )b h (D k )\D k ) = l^exp|-^i}^^/^£>|; s ,(fi). 

•h%(Q, j)n* k (n, j). (35) 
In these expressions the spin-rotated and shifted overlap-functions are given by 

n lk (n, j) ee dct Ns X* k (n, j), (36) 

and 

n%(n, j) ee Yl \ X h,h-{n, j)}- 1 Xl k ip (n, j) , (37) 

h'e\D>) 

respectively. The determinant in eq. i|36fl has to be taken for the occupied N e x A^ e -dimensional part of the (2 N x 2 N)— 
matrix X %k 

Xab = 2li ^cur, a^aa, a«r'(^> j)^aa-',b (38) 



4v 
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and in eq. I|37|) the indices h and p run over all the occupied and unoccupied states of the HF-determinant \D k ), 
respectively. (AT lfc ) _1 denotes the inverse of the occupied part of the matrix (38) and 



S aa , aa >(Sl,j) = D l Jl,{tt) exp 



2tt 
N 



In eq. i?>?>\) the spin-rotated and shifted matrix elements of the Hamiltonian involve also 



where 



and 



h l *(n, j) 



eta, ■ycr' 



, 2tt 

— ztcos — a 
N 



(fi, j) 



.'2tT 

-2tcos — a 
N 



+ ~~M a+f3-~f-SP5-<T,f3-o-{ "l J) 



■ IS 



N 



^, w (n,i) = ^V,y(n, i)]- 1 ^,^ 



Furthermore, introducing the definitions 



(39) 



(40) 



(41) 



(42) 



(43) 



(44) 



for all the occupied states h in the Slater-determinant \D k ), and 



-hi 
w 7<r 



iP ,(n, i) = ]T[i - p fci (n, jl^w^^-^ffi, j)^<», P < 



(45) 



for all the unoccupied states p' belonging again to the transformation D k , one can write the spin-rotated and shifted 
energy function out of eq. lj35l) as 



(46) 



The above described variational procedure will be denoted by LMSPHF (from Linear Momentum and Spin Projected 
Hartree-Fock) in the following. 

Finally, it should be stressed that this procedure can be used for any number of electrons N e = 1, . . . ,2N. As we 
already mentioned, for odd N e the spin quantum numbers 5 and £ are half-integer numbers. 



C. Spectral functions 

Let us now assume that we have obtained the ground state for an even number of electrons N e on a lattice with N 
sites solving the variational equations out of the sec. Ill Bl As we shall see in sec. Hill this ground state has always total 
spin 5 = but not necessarily linear momentum quantum number £ = 0. For instance, in the case of half-filling 
and N/2 even the ground state is obtained for £o = N/2. The underlying HF-transformation will be denoted by D 1 
and the total energy by Eq. For 5 = 0, there exists only a single coefficient fo, which is uniquely determined by the 
normalization 



(47) 
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The projected ground state of the iV e -electron system can thus be written as 

ID 1 ; N e £ S = E = 0) = P ° i0 C(6) ID 1 ) < 1/2 . 
We shall now approximate the (N e — l)-electron system (characterized by S = 1/2 and by 

m 1/2 

|ft;JV e -U_ lJ Sr=l/2cr) £ ^C^-iM^PV^, (49) 

i = l h=l ct'=-1/2 

where ID 1 ) refers to the HF-determinant for the iV e -electron ground state solution, while \D l ) for i = 2, . . . ,m refer 
to the determinants obtained for the lowest m — 1 excited states (which may correspond to different spin- and linear 
momentum quantum numbers than the ground state one's). The mixing coefficients / are then obtained by solving 
a generalized eigenvalue problem similar to eq. I|18|l 

(H - E~ h N)f = (50) 

with the constraint 

pNf = l 2m . Ne . (51) 
The overlap- and Hamiltonian-matrices are given by 

1 N 

N ihrT>kh , a , = - ^exp{-ifC-ii}8^/^D 1 ^,(0)- 
iV j= i 

■n ik (n,j)ni k h ,(£l,j) (52) 



and 

1 N 

H iha , kh , a , ee - Ew{-iwZ-ij}&J<K lDl J,l>W- 

■ tf k (n,j)hf jh ,(n,j), (53) 

where the short-hand notations in the integrands are defined by 

n% h ,(il,j) ee [X% th (il, j)]- 1 (54) 

and 

hf th ,(n, j) ee [x% th (ii, i)]- 1 ^^, j) £ i)fSr,^(n, i) . (55) 

respectively. Here 

*£, lh (n, j) ee J2 s-ywfa j)D k ' a „jAx^. h {^ j)}- 1 , (56) 



and all the other functions are defined in sec. Ill Bl For each possible linear momentum quantum number £_i the 
eqs. (l5f)l) and {51]) yield 2m • N e solutions h with 5=1/2 and energies E^. 

Similarly, the N e + 1-electron system (characterized by S = 1/2 and £+i ) will be approximated by 

m 2N 1/2 

\p ; N e + 1 e+i S = 1/2 a) = ]T E E ^X^+O^D^D*) , (57) 

i=l p=AT e + l <r' = -l/2 

where again |D : ) refers to the HF-determinant for the 7V e -electron ground state solution, while the |D l ) for i = 
2, . . . , to are taken from the lowest to — 1 excited states. The mixing coefficients g are obtained by solving a generalized 
eigenvalue problem 

(H - EpN)g = (58) 
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with the constraint 

The overlap- and Hamiltonian matrices are given by 



(59) 



■n ik (n,j)nf tP ,(n,j) 



(60) 



and 



1 n 

- Eexp{-i^ +1 i}^/dnr»^(fi). 



f^ipa, kp' a' 

iv J = l 

■ n ik {Q,j)h P k p ,(Q,j), 

where the short-hand notations in the integrands are defined by 



T? p k pl {n, j) = x; k p ,(n, j) - ]T x ik h [x ik h ,{n, j)\- l x%, p ,{n, j) 

h,h> 



and 



h p k pl (Q, j) ee n p k pl (Q, j)V k (n, j) + i)f 



t, 7 ^(n,i)^, P '(^i)> 



respectively. Furthermore, 



j) EE ^ D^ ip [l - j)W 

0,a' 



(61) 



(62) 



(63) 



(64) 



and all the other functions are again defined in sec. Ill Bl For each possible linear momentum quantum number £+i 
the eqs. I|58|) and il5H|l yield 2m • (2N — N e ) solutions p with 5 = 1/2 and energies E p . 

Now the hole-spectral functions can be calculated. The spectral function for hole states, S^k, efj, is defined by 



^S- h (k, e h ) ee (h;N e -lZ-iS=l/2\\% o -e_ 1 \\D 1 ;N e t S = 0) 

m N e 1/2 

= - ^o 1/2 EE E C'C- 

t=l fc=l a=-l/2 

JV , ^ , 1/2 



.7=1 k J <r'=-l/2" / 



ff' = -l/2 

iV e 



A'=l 



,(n). 
n a (n, j), 



(65) 



where A: = (27r/iV)£_i and er = Eq — E^. Similarly, the particle spectral function can be obtained from 



y/Sp(k, e p ) ee (p;iV e + U+i5=l/2||4 +i _ Co || J D 1 ;7V e eo5 = 0) 



m 27V 1/2 



= -^ 1/2 E E E 
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1/2 f+i* 



i=l p=Ar c + i cr = _i/2 



i=l *• ' a' = -l/2 J 

2N 



(66) 
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where now k = (2ir/N)£ t+ i and ep = Ep — Eq, 
The sum 



2m- N e 

J2 S- h (k, e h ) = n(k = {2-k/N)^) (67) 

h=l 

gives the occupation number of the basis state (jHJ in the ground state of the 7V e -electron system. Here instead of 
using £_i = 0, ... ,7V — 1, we use the equivalent representation £_i = —N/2 + 1, . . . , N/2. 

Furthermore, for plotting the spectral functions versus the linear momentum k and excitation energy uj = (or = ep) 
it is useful to introduce some artificial width for each state in order to obtain continuous functions of the excitation 
energy. For this purpose we use a Lorentzian shape with fixed width of 0.05 t for each hole (or particle) state. This 
also simplifies the representation of the density of states 

N/2 

N{u>) = ]T [S (h) ((2n/N)a, to) + S (p) ((2ir/N)a, w)] . (68) 

a=-N/2+l 

where the indices h and p have been put here in parentheses, since they are absorbed now in the continuous variable 

Ul. 

In order to determine the hole spectral functions from an exact (Lanczos) solution, we start with the exact ground- 
state of the system with N e electrons, |7V e )o, which is represented by a linear combination of the basic configurations 
and apply the annihilation for an electron with spin projection a and momentum a = £o — £-i 

Ca,aWe)o = Va\N e - 1,£_ 1; S = l/2) X . (69) 

The resulting state \N e — 1, S = 1/2) x is a state with N e — 1 electrons and well defined quantum numbers for the 
momentum (£-i) and spin. The constant rj a has been introduced to normalize the state. This state is in general not an 
eigenstate of the Hamiltonian, but it can be used as a starting point to generate further states \N e — 1, S = 1/2) i 
with the same quantum numbers by means of the Lanczos method (T3|. A diagonalisation of the Hamiltonian in the 
space which is generated by these basis states \N e — 1, S = 1/2) i for i = 1 . . . v will lead to v eigenvalues e l a and 
eigenstates \e l a ), which converge with increasing v to the exact solutions for the system with N e — 1 electrons. The 
spectral function for the hole states can now be expressed as 

V 

S h ((2n/N)a,Lu)= V lY / \(£«\Ne-l,Z-i,S= 1/2)^5(^-6^ • (70) 
i=i 

This means that the occupation probability n((2-n: / N)a) is determined by the square of the normalisation constant r] a 
defined in eq. l(69"jl and the overlap of the states \e l a ) with the starting vector for the Lanczos iteration | N e — 1, S = 
1/2)! defines the spectral strength at the energy oj = e l a . The result for the spectral function converges very rapidly 
with the number of iterations v. A corresponding procedure starting with a state 

clMo, (71) 
can be used to determine the spectral function for the particle states. 



III. RESULTS AND DISCUSSION 



First, we checked the LMSPHF-approach discussed in sec. Ill Bl for half-filled lattices associated to N — 2, N = 4 
and N = 6 sites, respectively. Simple counting of the number of variational variables (see tab. H} shows that for 
these lattices the approach is exact. Each state can be represented by a single spin- and linear momentum-projected 
HF-configuration for any value of U/t. This has been tested explicitely for the 36 different states of the half-filled 
N = 4 lattice as well as for the 5 lowest states of each possible linear momentum quantum number with spins 5 = 0, 
S = 1 and S — 2 in the half-filled TV = 6 grid. As an illustrative example we present in Fig. the energy spectrum 
for the half-filled N = 4 lattice as obtained for U/t — 4. Each of these states can be represented by a single projected 
LMSPHF-determinant. The figure also shows nicely the degeneracy of the states with £ = 1 and £ = 3 (or -1) 
discussed in sec. Ill Al As expected, the ground states have all total spin 5 = 0, however, while for the lattices with 
N = 2 and N = 6 the linear momentum quantum number £ = is obtained, the TV = 4 ground state has £ = TV/2 = 2. 
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TV 



4 



6 12 14 



30 



Hilbert restricted space 
variationnal parameters 



~ 10 



~ 70 - 2 x 10 4 ~ 2 x 10 1 



,5 



~ 10 



32 



72 



288 



392 



1800 



TABLE I: Comparison of the dimension of the Hilbert restricted space (taking account of the symmetries) for a given half-filled 
lattice with N sites with the number of parameters involved in the LMSPHF-approach. 

This feature persists, if half-filled lattices with larger TV are considered : For TV/2 being an odd integer the ground 
state has 5 — and £ = 0, while for N/2 being an even integer S = and £ = N/2 is obtained. 

This effect can be understood by a simple spin-correlation. Assume for the moment small interaction strength U. 
Then, if N/2 is odd, the configuration with the lowest energy is obviously the determinant in which the single particle 
states (jHJ with a — 0, ... , ±(TV — 2)/4 are all filled, each by two electrons with opposite spin-directions. The resulting 
total spin is then 5 = and the linear momentum quantum number £ = 0. For N/2 being even, however, the situation 
is different. Here, the last two electrons have to be distributed over the two degenerate states with a = ±iV/4. 
It can be shown that for two electrons with opposite spin directions occupying two different basis orbits <PJ always the 
5=1 component is energetically favored with respect to the 5 = component. The spin-correlation explains 
why half-filled lattices with N/2 being odd has the lowest excited state associated with total spin 5 = 1 and 
linear momentum quantum number £ = N/2, Assuming again small interaction strength U, the energetically lowest 
configurations are obtained by promoting one of the four electrons from the last occupied orbits (a = ±(TV — 2)/4) to 
the first unoccupied orbits (a — ±(N + 2)/4) in the above mentioned (5 = 0, £ = 0)-ground state configuration for 
odd N/2. There are 8 degenerate determinants of this type, 4 of them with £ = N/2, 2 with £ = 1 and 2 with £ = — 1 
(or, equivalent £ = N — 1). The expectation value of the Hamiltonian for each of these determinants will be denoted 
by b. The £ = ±1 configurations yield two degenerate S = 1 states at energy b — U/N and two degenerate 5 = 
states at energy b + U/N. The four £ = N/2 determinants, however, can be coupled to two 5=1 configurations both 
with energy b — U/N and an interaction of —U/N between them and to two 5 = configurations both with energy 
b + U /N and an interaction of U /N in between. Thus we get here one 5=1 solution with energy b — 2U /N, an 
5=1 and an 5 = solution both having energy b, and one 5 = state with energy b + 2U/N. Consequently, here 
the lowest (5 = 1, £ = TV/2)-solution will become the first excited state. Again the relative splitting in between the 
various states will even increase with increasing interaction. Also this expectation is confirmed by the results obtained 
for the half-filled TV = 6-, N = 14- and TV = 30-lattices in the present work. 



Let us now consider the half-filled TV = 12 lattice. Here in the usual Lanczos approach (all £ = determinants) 
already 853 776 configurations have to be treated, and even using all the symmetries, for 5 = and the 12 possible 
£-values still between 18 840 and 18 916, for 5 = 1 and the possible £-values between 31 833 and 31 872 configurations 
have to be accounted for. Fig.[2displays the energies of the first excited states for various methods using an interaction 
strength of U/t = 4. In the rightmost column (EXACT) the results of the Lanczos diagonalisation for the lowest 4 
states are presented, the other columns refer to different variational approaches. The leftmost one (HF) gives the 
Hartree-Fock result obtained by using the determinant (12) as test wave function and not caring about spin- and 
linear momentum. The next two columns have been obtained by HF calculations with linear momentum projection 
before the variation only (LMPHF) but still not restoring the total spin and with spin-projection before the variation 
only (SPHF) but not restoring the total linear momentum. The next column (LMS Z PH) results from variational 
calculations with projection on good linear momentum and z-projection of the spin S z = 0. Finally, the second 
but last column (LMSPHF) presents the results obtained for the lowest 5 states with linear momentum- and full 
spin-projection before the variation as described in sec. Ill Bl As can be seen, the LMSPHF-approach reproduces the 
energies of the exact solutions not only for the ground but also for the lowest excited states very well. The deviations 
vary only between 0.16 and 0.58 percent. This is remarkable since the number of variational parameters (288 for 
5 = in this case) is significantly smaller than the dimension of the space with good symmetries, which is 18840 in 
this example. 

As expected from the above arguments based on spin-correlation the ground state has 5 = and £ = TV/2 = 6, while 
for the quantum numbers of the first excited state 5=1 and £ = are obtained. It is furthermore obvious, that 
both, linear momentum and spin have to be restored simultaneously. All results obtained by performing none or only 
part of the corresponding projections fail to reproduce the exact spectrum. 

Fig. presents the energies of lowest 5 = 0, 5 = 1 and 5 = 2 states obtained with the LMSPHF-approach for the 
various linear momentum quantum numbers £. Note, that the £ = 1-results are degenerate with those for £ = 11, 
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Half-filled TV=12 lattice 
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those for £ = 2 with those for £ = 10, etc., so that only the results for £ = 0, ... 6 are shown. Except for the 
(5 = 0, £ = ./V/ 2) -ground state, the lowest states for all other linear momentum quantum numbers have total spin 
5 = 1. This supports again the importance of the above discussed spin-correlation. 

In fig. El we represent the occupation numbers of the various basis orbitals (jSJl in the LMSPHF-ground state obtained 
via eq. IJ67H and we compare it with those resulting from the one-body density of the exact ground state and with 
those expected for a non-interacting Fermi-gas. Again, excellent agreement of the LMSPHF-results with the exact 
solution is obtained. It should be stressed, that the occupation numbers obtained via eq. ll67T) are not affected by the 
number of determinants m used for the calculation of the hole-spectral functions. Identical numbers are obtained, no 
matter whether only the ground state determinant \D l ) (m = 1) or, e.g., all the m = 5 determinants corresponding 
to the five lowest LMSPHF-solutions Fig. [21 are included in the calculation. 

In fig. we compare the occupation numbers of the basis orbitals in the LMSPHF-(5 = 0, £ = 6)-ground states 
obtained for various strength parameters U ft of the interaction. Since in all cases the LMSPHF- and exact occupation 
numbers cannot be distinguished, we have plotted only the former. As expected, the correlations (i.e, the deviation 
from the Fermi-gas values) grow with increasing interaction. Already at U/t = 64 the result looks rather similar to 
the equal distribution of the anti-ferromagnetic limit expected for U/t — ► oo. 

In fig. El we represent the hole- (eq. I|65|) ) and particle- (eq. (jSEl) spectral functions versus the excitation energy 
to (in units of t). Here, to = for the hole- and to = ep for the particle-states, respectively. We compared the 
results obtained by using only one determinant (m = 1, upper half of the figure) with those resulting from using 
the HF-transformations obtained for all the m = 5 lowest LMSPHF-states presented in Fig. [51 As expected for the 
half-filled lattice, particle- and hole-states are nicely symmetric around the Fermi-energy top = U/t/2. They are 
separated by the so-called "Hubbard-gap" of about U/t/2. Using m = 5 determinants one obtains 120 hole and 120 
particle states for each possible linear momentum quantum number, while we obtain only 24 states for m = 1 for 
each £±i-value. Consequently, the strength is more spread for the m — 5-approximation with respect to the more 
restricted m = 1-approximation. 

Finally, in Fig. we compare the results for the strength defined by eq. ilfifil and obtained with m = 5 determinants 
for 5 — 1/2 and £ = 0, ... ,3 with those resulting from the exact calculation as described in l|69"j) - (ITDl . As can be 
seen, the agreement for both the excitation energies as well as for the splitting of the strength is again very good. 

B. Doped lattice 

Let us now consider a system of only 10 electrons in the N = 12-lattice, again using U/t = 4 as interaction strength. 
N e /2 is odd, then all the basis orbits (EJ with a = 0, ±1 and ±2 can be filled, each with two electrons with opposite 
spin directions. Thus, from spin-correlation argument, one expects an (5 = 0, £ = 0)-ground state. The lowest 
excited states can then be obtained by promoting one of the four electrons from the a = ±2-orbitals to the a = ±3 
orbits. There are again 8 degenerate determinants of this type : 2 with £=1,2 with £ = — 1, 2 with £ = 5 and 2 
with £ = — 5. Calling C the expectation value of the Hamiltonian © within each of these determinants, one obtains 
four degenerate 5 = 1 states at energy C — U/N and four degenerate 5 = states at energy C + U/N, respectively. 
The corresponding linear momentum quantum numbers are in both cases £ = 1,-1,5 and —5. Thus one expects the 
lowest excited states to have 5=1 and £ = ±1, ±5. 

Fig. [SI presents the energies of lowest 5 = 0, 5 = 1 and 5 = 2 states of the 10-electron system on the N = 12-lattice 
obtained with the LMSPHF-approach for the various linear momentum quantum numbers £ = 0, ... 6. As in Fig. 
the spectra for £ = —1, . . . , —5, which are degenerate with those for £ = 1, ... 5 have not been plotted. As expected, 
we obtain an (5 = 0, £ = 0)-ground state, and the lowest four excited states have all 5 = 1 and linear momentum 
quantum numbers of £ = ±1 and £ = ±5, respectively. 

Obviously, the nice particle-hole symmetry in the spectral functions of the half-filled N = 12-lattice is destroyed, if 
only 10 electrons in this lattice are considered. Using the simplest (m = 1) approximation for the calculation of the 
5 = 1/2-wave functions for the 9- and 11-electron systems, one obtains here for each possible linear momentum 20 
one-hole- and 28 one-particle-states. This is reflected in the hole- and particle-spectral functions, which are shown 
in Fig. Has a function of the excitation energy. Particle- and hole-strengths are now distributed asymmetrically 
around the Fermi-energy and furthermore, the Hubbard-gap obtained for the half-filled lattice has vanished. 

C. N=14 lattice 

We shall now consider the half-filled N = 14-lattice, again using U/t = 4. Here, in the usual Lanczos-approach (all 
determinants with S = 0) already 11 778 624 determinants have to be treated. Consequently, for this example we 
could obtain only the exact ground and first exited state within about one week of computer time. Even if all the 
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symmetries would be used, for 5 = and the 14 possible linear momentum quantum numbers still between 197 099 
and 197 276, for 5 = 1 and the possible ^-values between 357 770 and 357 945 configurations have to be accounted 
for. 

In fig. 1101 we display the results for the lowest states obtained with different variational methods and compare them 
with the exact results. The nomenclature is the same as in Fig. |21 Again it is seen, that all results obtained by 
performing none or only part of the symmetry-projections fail to reproduce the exact spectrum. Only if linear- 
momentum and full spin-projection are both performed simultaneously before the variation as done in the LMSPHF- 
approach the exact data can be reproduced. The deviation of the LMSPHF-energies from the latter amount here to 
0.38 percent for the ground and 0.93 percent for the first excited state, respectively. Thus (as expected because of 
the considerably larger dimensions) they are larger than in the half-filled TV = 12-example. It should be stressed, 
however, that in the LMSPHF-approach each of these states are represented by only a single symmetry-projected 
configuration. Correlating these solutions by additional determinants could, obviously, still improve the agreement. 
As expected for a half-filled lattice with TV/2 being odd, the ground state has total spin 5 = and linear momentum 
quantum number £ = 0, while for the first excited state 5 = 1 and £ = N/2 = 7 is obtained. 

In fig.Elwe represent the energies of the lowest 5 = 0, 5 = 1 and 5 = 2 states obtained with the LMSPHF-approach 
for the various linear momentum quantum numbers £. Again, the spectra for £ = TV/2 + 1, . . . , (N — 1) are degenerate 
with those obtained for £ = 1,.. . N/2 — 1, respectively, and are hence not displayed. Except for the (5 = 0, £ = 0)- 
ground states and the £ = 2 case, where the lowest 5 = and 5 = 1 states are almost degenerate, the lowest states 
for all the other linear momentum quantum numbers have always total spin 5 = 1. This supports the above discussed 
spin-correlation favoring the 5 = 1-excitations. 

Also for the half-filled TV = 14-lattice, the LMSPHF- and exact occupation numbers of the ground state can hardly 
be distinguished and are not shown here. Instead, we present in Fig. El the hole- and particle-spectral functions 
obtained with the m = 5 HF-transformations resulting from the five lowest LMSPHF-solutions out of Fig. El In 
order to obtain continuous functions of the excitation energy w = e^ (or u> = eg), each state has been broadened with 
a Lorentzian of constant width 0.05 t. As can be seen, the particle- and hole-strengths are symmetrically distributed 
around the Fermi-energy (again ujf = U/t/2) and separated by a Hubbard-gap of about the same size (U/t/2) as 
obtained for the half-filled TV = 12-lattice. 



D. Large TV=30 lattice 

Finally, we shall report LMSPHF-results for the half-filled TV = 30-lattice, again for interaction strength U/t = 4. 
Here in the usual Lanczos-approach 2.9 • 10 15 states would have to be included, which is obviously impossible. Even 
coupling the configurations to good total spin and linear momentum quantum number, still for the various ^-values 
dimensions of the order of 9.7 • 10 13 for 5 = and even 2.2 • 10 14 for 5 = 1 would have to be treated. In the 
LMSPHF-approach, however, for each 5 = 0-state "only" functions of 1800 real variables have to be minimized. 
As in Figs. Eland El we present the results for the energies of the lowest states as obtained with various variational 
methods for the half-filled TV = 30-lattice in Fig. El Obviously, there is no exact result available to compare with. 
According to the spin coupling arguments, since TV/2 is odd the ground state has 5 = and linear momentum 
quantum number £ = while the first excited state is obtained for 5 = 1 and £ = TV/2 = 15. Again it is seen that the 
simultaneous restoration of linear momentum and total spin before the variation as done in the LMSPHF-approach 
is essential. All other approximations produce results, which do not come close to the LMSPHF-energies. 
The occupation numbers for the different basis states (|SJl in the LMSPHF-ground state are presented in comparison 
with the values expected for a non-interacting Fermi-gas in Fig. El Like for the half-filled TV = 12 lattice (see Fig. El, 
strong effects of the correlations are observed. 

The density of states N(cu) (eq. jSHJ) for the half-filled TV = 30-lattice is shown in Fig. El Here, we have performed a 
simple 1-determinant calculation (m = 1) as explained in sec. Ill CI Again, the states were broadened by a Lorentzian 
with a width of 0.05 t. As expected for half-filling, particle- and hole- strengths are distributed symmetrically. The 
Fermi-energy is again uop = U/t/2 and the width of the Hubbard-gap again of about the same size as obtained for 
the half-filled TV = 12- and TV = 14-lattices (see figs. El and El- • 

Finally, in Fig. El the particle- and hole-spectral functions are presented as functions of linear momentum k and 
excitation energy w. Again, Lorentz-shape and a constant width of 0.05 t has been used. Obviously, the strength 
of each of these spectral functions would still be redistributed, if instead of only one, several determinants would be 
taken into account for the calculation. Since anyhow, there are no "exact" (or experimental) data to compare with, 
this generalization has not been done for the TV = 30-lattice in the present investigation. 
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IV. CONCLUSIONS 

We have devised a variational approach for the approximate solution of the 1-dimensional Hubbard-model with 
nearest neighbor hopping and periodic boundary conditions. For the ground state we start with a Hartree-Fock type 
transformation mixing all the quantum numbers of the single electron basis states. 
The results may be summarized as follows : 

For many-electron systems on small lattices, where the number of configurations to be treated for each pair of good 
spin- and linear momentum-quantum number is smaller than the number of real variational degrees of freedom 
accounted for in the variational calculation, the LMSPHF-approach is obviously exact. However, even in the half- 
filled N = 12- and N = 14-lattices, where for each pair of quantum numbers (S, £) the number of configurations is by 
2 or even 3 orders of magnitude larger than the number of variational degrees of freedom, the LMSPHF-approach does 
reproduce the exact energies of the ground and lowest excited states with an accuracy of better than 99 percent, the 
occupation numbers of the ground and lowest excited states even better, and yields (at least if several determinants 
for the calculation of the one-hole and one-particle spectra are included) even for the spectral functions very good 
agreement with the exact results. This gives some confidence, that the LMSPHF-approach can be considered as a 
very good truncation scheme even for lattices which are far too large to allow for exact diagonalisation. Even where 
complete diagonalisation via the Lanczos method is still numerically possible, the LMSPHF-approach is much faster. 
Spin- and linear momentum- quantum numbers for the ground and the lowest excited state can reliably be predicted 
making use of the fact that for two electrons with opposite spin-directions in two different momentum space basis 
orbits always the S = 1-configuration is energetically favored with respect to the S = 0-one. Thus, for half-filled 
lattices with N/2 being odd the ground state has S = and linear momentum £ = 0, while the lowest excited state 
has 5=1 and £ = N/2. For half-filled lattices with N/2 even on the other hand one obtains S = and £ = N/2 
for the ground and S = 1 and £ = for the first excited state. The same spin-correlation can be used to predict the 
quantum numbers of the lowest states away from half-filling and is also supported by the observation that for almost 
all linear momentum quantum numbers (except for that of the ground state) the lowest state has always total spin 
S = l. 

One of the nicest features of the LMSPHF-method is that it still can be improved rather easily. So, e.g., instead of 
HF-determinants generalized Slater-determinants of the HF-Bogoliubov (HFB) type can be used as basic building 
blocks. This would increase the number of real variational degrees of freedom (for the example of an S = 0-state) 
from 2N e (2N — 2N e ) in the LMSPHF-approximation to 2N(2N — 1) for any number of electrons. The price one has 
to pay for this, is an additional integration due to the (then necessary) projection onto good total electron number 
in the calculation of all the matrix-elements. Because of the simple form of the interaction, however, this should 
not cause any serious problems. The resulting linear momentum-, spin- and number projected HFB-approach would 
consequently account for considerably more correlations in each single configuration than the LMSPHF-method. 
Furthermore, instead of using essentially only one symmetry-projected configuration for each state, correlating 
symmetry-projected configurations can be added and the underlying HF- (or HFB-) transformations again be de- 
termined by variational calculations. Since the energy gain obtained for a particular state under consideration due 
to the last added configuration is by construction always smaller than the energy gain due to the last configuration 
added before, this procedure (which already has been applied successfully to the nuclear many-body problem jjj]) can 
also give a hint on the quality of the solution even in lattices, which are too large to allow for an exact solution. 
Last but not least, the procedure discussed in the present work can be extended easily to the 2-dimensional Hubbard- 
model, which is supposed to be of larger physical relevance than its 1-dimensional simplification. Then, obviously, the 
total linear momentum becomes a 2-dimensional vector and one has to project on each of its components separately. 
Thus again, with respect to the LMSPHF-approach discussed above only one additional integration is needed. 
This leaves ample space for future investigations and we are quite confident that the variational approaches originally 
devised for the nuclear many-body problem will turn out to be rather useful even for more complicated lattice Fermion 
models than the simple 1-dimensional version investigated in the present work. 
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FIG. 1: Energy spectrum of the half-filled N = 4-lattice as obtained with interaction strength U/t = 4 either by complete 
diagonalization or by the LMSPHF-approach out of sec. Ill Bl In the latter case each of the states is represented by a single 
linear momentum- and spin-projected Slater-determinant. 
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FIG. 2: Energies for some states of the half-filled iV = 12-lattice as obtained using the interaction strength U/t = 4 by various 
approaches. The different columns refer to a simple unprojected HF-calculation (HF), to HF with only linear momentum- 
(LMPHF), or only spin-projection (SPHF) before the variation and to HF with restoration of the total linear momentum 
and only the 3-component of the total spin (LMS Z PHF). Finally, the energies of the lowest few states obtained with linear 
momentum- and full spin-projection before the variation as described in sec lll Bl fLMSPHF) are compared with those resulting 
from a complete diagonalization (EXACT). 
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FIG. 3: Again for the half-filled iV = 12-lattice and U/t = 4, the energy spectra for the lowest S = 0, S = 1 and S = 2 states 
as obtained with the LMSPHF-method for the various linear momentum quantum numbers £ are displayed. Because of the 
degeneracy of the £ = 12 — i with the £ = i spectra for i = 1, . . . , 5, only the spectra for £ = 0, . . . 6 are presented. 




FIG. 4: For the half-filled iV = 12-lattice and U/t — 4 the occupation numbers (eq. (jlHJl) of the various basis states in the 
(S = 0, £ = 6)-LMSPHF-ground state are compared to those obtained by a complete diagonalization (EXACT) and to those 
expected for a non-interacting Fermi-gas. Identical occupation numbers are obtained for using only the lowest (m = 1) or the 
lowest five (m = 5) determinants in the calculation of the hole-spectral functions. 
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FIG. 5: Same as in Fig. 0] but for various interaction strengths U/t. Since in all cases the exact and LMSPHF-results cannot 
be distinguished, only the latter are displayed. 
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FIG. 6: Again for the half-filled N = 12-lattice and U/t = 4, the hole 
plotted for all possible values of £±i as functions of the energy uj ( 

has been obtained by using only the ground state determinant ID 1 ) in the calculation for the spectral functions (m = 1). Using 
instead all the m — 5 determinants corresponding to the lowest five states out of Fig. [2] in these calculations, one obtains the 
results displayed in the lower part of the figure. 
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FIG. 7: Again for the half-filled N = 12-lattice and U/t = 4, the hole spectral functions for £ = 0, . . . , 3 as obtained 
with m = 5 determinants are plotted versus the excitation energy and compared with the "exact" results computed with the 
Lanczos-approach as discussed in the text. 
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FIG. 8: Same as in Fig. El but for only 10 electrons in the N = 12-lattice. 
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FIG. 9: Same as in Fig.HjJ but for only 10 electrons in the A*' = 12-lattice. Here only the results obtained with one determinant 
(m = 1) are presented. 
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FIG. 10: Same as in Fig.EI but for the half-filled N = 14-lattice. 
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FIG. 11: Same as in Fig.H but for the half-filled N = 14-lattice. 
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FIG. 12: Again for the half-filled N = 14-lattice, the sum of the hole- and the particle-spectral function are displayed as 
functions of the excitation energy tJ (in units of t) and the linear momentum k in units of (2n/N). The spectral functions have 
been obtained here by using the m — 5 determinants corresponding to the five lowest LMSPHF-solutions in Fig. ED In order 
to obtain continuous functions for each state a Lorentz-shape with constant width of 0.05 t has been used. 
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FIG. 13: Same as in Fig. El but for the half-filled N = 30-lattice. A gain U/t = 4. Here no exact results are available to 
compare with. 
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FIG. 14: Same as in Fig. 0] but for the half-filled iV = 30-lattice. Again U/t = 4, and no exact results are available. 
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PIG. 15: For the half-filled iV = 30-lattice with U/t = 4 the total number of states N(u>) out of eq. I(i8ll is presented as function 
of the excitation energy. Here only m = 1 determinant was used to obtain the one-hole- and one-particle states. Each of these 
states has been artificially broadened to a Lorentz-shape with constant width of 0.05 t. 




FIG. 16: (Color online) Same as in Fig. [121 but for the half-filled N = 30-lattice, and using only m = 1 determinant in the 
calculation of the spectral functions. 



